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Charm production is studied in deep-inelastic electron-photon scattering using 
OPAL data at e + e - centre-of-mass energies from 183 to 209 GeV. Charm quarks 
are identified by exclusive reconstruction of D* mesons. The cross-section of D* 
production is measured in a restricted kinematic region, and then extrapolated to 
the total charm production cross-section and the charm structure function of the 
photon. For x > 0.1 the measurement is well described by Monte Carlo models 
and perturbative QCD calculations but for x < 0.1 the predictions are lower than 
the data both in the directly measured region and after the extrapolation. 


1 Introduction 

The charm component of the photon structure function, c , has been mea¬ 
sured by OPAL at LEP2 by applying the well established method of exclusive 
D* reconstruction to deep-inelastic electron-photon scattering events. The 
determination of F% c exploits the fact that the differentialmross-section as 
function of Q 2 and Bjorken x is proportional to F^ c {x, Q 2 ). Ill 

Due to the large scale established by their masses, the contribution to F% 
from charm quarks can be calculated in perturbative QCD, and predictions 
have been evaluated U at next-to-leading order (NLO) accuracy. F ^ c receives 
contributions from the point-like and hadron-like components of the photon 
structure, with the hadron-like component dominating at very low values of 
x and the point-like part accounting for most of F^ c for x > 0.1. 

The preliminary results presented here a extend the earlier measurement u 
of E7 c using basically the same analysis strategy. It is based on 654.1 pb _1 
of data for e + e _ centre-of-mass energies from 183 to 209 GeV, recorded by 
the OPAL experiment in the years 1997-2000. 


* Presented at the International Conference on the Structure and Interactions of the 
Photon, including the 14th International Workshop on Photon-Photon Collisions, Ascona, 
Switzerland, 2-7 September, 2001. 
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Figure 1. Distribution of the difference between the D* and D° candidate masses. The data 
are shown as points with statistical errors, while the histogram represents the combinatorial 
background estimated using events with wrong-charge combinations for the decay products 
of the D* mesons. The curve is the result of the fit to the data. 


2 Data selection 

The most important cuts used in the selection of deep-inelastic electron- 
photon scattering events containing a D* are summarised below. 

1. An electron candidate must be present with an energy E ta g > 0.5-Eb and 
a polar angle in the ranges 33 < 9 tag < 55 mrad (SW) or 60 < 0 tag < 
120 mrad (FD), corresponding to 5 < Q 2 < 100 GeV 2 . 

2. Double-tag events are eliminated by requiring that the sum of all energies 
in the SW and FD detectors opposite to the tag are below 0.25^ 

3. An exclusively reconstructed D* candidate must be present with a trans¬ 
verse momentum p® > 1(3) GeV for SW(FD)-tagged events and a pseu¬ 
dorapidity |? 7 D | < 1.5. The D* meson must decay into D°7 t with the D" 
decaying into the charged particles K7r or K7T7r7r. 

Figure [I] shows the difference between the D* and D° candidate masses for 
both decay channels combined. A clear peak is observed around 0.145 GeV, 
the mass difference between the D* and the D° mesons. An unbinned max¬ 
imum likelihood fit to this distribution gives 60.3 ± 10.3 signal events above 
the combinatorial background from deep-inelastic electron-photon scattering 
events e + e - —> e + e _ qq with q=uds. The expected background from all other 
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Figure 2. Data distributions compared to the HERWIG and Vermaseren predictions. For 
HERWIG several predictions are shown: the full prediction, the point-like component alone 
(HW PL), and a superposition of the HERWIG point-like prediction together with a scaled 
hadron-like prediction, denoted by HW scaled. 


processes that potentially contain D* mesons in the final state is found to be 
negligible using Monte Carlo simulations. 

Figure || shows the distributions of two global event quantities, Q 2 and 
W v i s , and two variables related to the kinematics of the D* candidates, p® and 
|77 D * |. The data are compared to the absolute predictions of the HERWIG6.I □ 
and Vermaseren □ leading order Monte Carlo programs. 

To get a better description of the data, the hadron-like component of 
the HERWIG prediction has been fitted to the = 2p® /W v ; s distribution 
shown in Figure while keeping the point-like part fixed, resulting in a scale 
factor of 6.6 ± 2.7. There are several possible sources for this difference. The 
NLO prediction itself has a significant uncertainty due to variations of the 
charm quark mass and the renormalisation and factorisation scales, the gluon 
distribution of the photon has large experimental errors, and uncertainties in 
the shape and the modelling of the p ^ distribution can change the efficiency 
for the selected events. 
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Figure 3. The measured x Jjj! distribution compared to the predictions described in the 
caption of Figure pi 


This scaled HERWIG prediction, also shown in Figure ||, is used to es¬ 
timate the signal selection efficiency. The difference between the results ob¬ 
tained with the scaled and the original HERWIG models is taken into account 
as a systematic uncertainty. 

3 Results 

Table [I] summarises the cross-section of D* production in deep-inelastic 
electron-photon scattering measured in the kinematic region defined by the 
event selection. The total cross-section for cc production in deep-inelastic 
electron-photon scattering, shown in Figure |]a), is the result of an extrapo¬ 
lation to the whole kinematic region using the HERWIG scaled model. The 
value of the charm structure function of the photon, Ej c (x, (Q 2 )) / a em , av¬ 
eraged over the corresponding bin in x, shown in Figure f|b), is obtained 
using the ratio c (x, {Q 2 ))/a em _/a{e + e~ —► e + e“ccX) given by the NLO 
calculation, i 


Table 1. The cross-section cr D measured in the restricted region, compared to Monte Carlo 
predictions. The numbers in parentheses refer to the point-like and hadron like components. 



0.0014 < x < 0.1 

0.1 < x < 0.87 

OPAL 
HERWIG 
HW scaled 
Vermaseren 

4.7 ±1.3 ±0.9 
1.02 (0.65 ± 0.37) 
3.10 (0.65 ± 2.45) 
0.84 

3.0 ±0.9 ±0.4 
2.05 (2.02 ± 0.03) 
2.23 (2.02 ± 0.21) 
2.81 
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Figure 4. Preliminary OPAL results for a) ir(e + e _ —> e+e _ cc ATkiwith 5 < Q 2 < 100 GeV 2 
and b) c (x, (Q 2 ))/a em . The band for the NLO calculation El indicates the theoretical 
uncertainties. 


All models and predictions shown in Figure |] are in good agreement with 
the measurement for x > 0.1, where the purely perturbative point-like process 
is dominant and both the experimental and the theoretical uncertainties are 
moderate. 

On the other hand, for x < 0.1 the measurement lies more than two stan¬ 
dard deviations above the predictions, which contain an approximately equal 
contribution from the point-like and hadron-like processes. This difference is 
also observed in the comparisons to the Monte Carlo models in the directly 
measured kinematic region, where the extrapolation uncertainties are avoided. 
Unfortunately, no theoretical calculations are available for comparison in the 
restricted phase-space. 
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